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We report quantitative analysis of positive magnetoresistance (PMR) for unidirectional-lateral- 
superlattice samples with relatively small periods (a=92-184 nm) and modulation amplitudes 
(Vb=0. 015-0. 25 meV). By comparing observed PMR's with ones calculated using experimentally 
obtained mobilities, quantum mobilities, and Vb's, it is shown that contribution from streaming 
orbits (SO) accounts for only small fraction of the total PMR. For small Vo, the limiting magnetic 
field B c of SO can be identified as an inflection point of the magnetoresistance trace. The major 
part of PMR is ascribed to drift velocity arising from incompleted cyclotron orbits obstructed by 
scatterings. 
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I. INTRODUCTION 

Large mean free path (L^>1 /im) of GaAs/AlGaAs- 
based two-dimensional electron gas (2DEG) and modern 
nano-fabrication technologies have enabled us to design 
and fabricate 2DEG samples artificially modulated with 
length scales much smaller than L. The samples have 
been extensively utilized for experimental investigations 
of novel physical phenomena that take place in the new 
artificial environments^ Unidirectional lateral superlat- 
tice (ULSL) represents a prototypical and probably the 
simplest example of such samples; there, a new length 
scale, the period a, and a new energy scale, the amplitude 
Vo , of the periodic potential modulation are introduced to 
2DEG. These artificial parameters give rise to a number 
of interesting phenomena through their interplay with 
parameters inherent in 2DEG, especially when subjected 
to a perpendicular magnetic field B. Magnetotrans- 
port reveals intriguing characteristics over the whole span 
of magnetic field, ranging from low field regime domi- 
nated by semiclassical motion of electrons through 
quantum Hall regime where several Landau levels are 
occupiedj£iL2i2ii& up to the highest field where only the 
lowest Landau level is partially occupie d 1 ^ 12 ' 13 ; in the 
last regime, semiclassical picture is restored with compos- 
ite Fermions (CFs) taking the place of electrons. Of these 
magnetotransport features, two observed in low fields, 
namely, positive magnetoresistance 4 (PMR) around zero 
magnetic field and commensurability oscillation 2 * 3 " (CO) 
originating from geometric resonance between the pe- 
riod a and the cyclotron radius R c =Tik-F / e\B\, where 
kp=\/2irn e represents the Fermi wave number with n e 
the electron density, have the longest history of being 
studied and are probably the best-known. 

The PMR has been ascribed to channeled orbit, or 
streaming orbit (SO), in which electrons travel along the 
direction parallel to the modulation (y-direction), being 
confined in a single valley of the periodic potential, 4 Elec- 
trons that happen to have the momentum perpendicular 
to the modulation (cc-direction) insufficient to overcome 
the potential hill constitute SO. In a magnetic field B, 



Lorentz force partially cancels the electric force deriv- 
ing from the confining potential. Therefore the number 
of SO's decreases with increasing B and finally disap- 
pears at the limiting field where Lorentz force balances 
with the maximum slope of the potential. The extinc- 
tion field B c depends on the amplitude and the shape of 
the potential modulation, and for sinusoidal modulation 
Vo cos(27rx/a), 



B r = 



2nm*Vn 



(1) 



where m* represents the effective mass of electrons. It 
follows then that Vq can be deduced from experimen- 
tal PMR provided that the line shape of the modula- 
tion is known, once B is determined from the analysis 
of the experimental trace. An alternative and more fa- 
miliar way to experimentally determine Vb is from the 
amplitude of CO. In the past, several groups compared 
Vo's deduced by the two different methods for the same 
samples ^i^iiSiii In all cases, Vb's deduced by PMR and 
by CO considerably disagree, with the former usually 
giving larger values. Part of the discrepancy may be at- 
tributable to underestimation of Vo by CO, resulting from 
disregarding the proper treatment of the decay of the 
CO amplitude by scatteringpiS* 1 ^ 2 ^ However, the most 
serious source of the disagreement appears to lie in the 
difficulty in identifying the position of B e from an ex- 
perimental PMR trace, which was taken, on a rather ad 
hoc basis, as either the peakfi 4 * 1 ^ 1 ^ or the position for 
steepest sloped It is therefore necessary to find out the 
rule to determine the exact position of B e . This is one of 
the purposes of the present paper. We will show below 
that B D can be identified, when Vb is small enough, as an 
inflection point at which the curvature of PMR changes 
from concave down to concave up. Another target of the 
present paper is the magnitude of PMR. The magnitude 
should also depend on Vq as well as on other parame- 
ters of ULSL samples. The subject has been treated in 
theories by both numericalSi* 2 ^ and analytical^ 3 , calcula- 
tions. However, analyses of experimental PMR is so far 
restricted to the qualitative level 4 " that the magnitude 
increases with Vq- To the knowledge of the present au- 



TABLE I: List of samples 



(a) 



No. 


a (nm) 


Hall-bar size (/^m 2 ) 


back gate 


1 


184 


64x37 


X 


2 


184 


64x37 


X 


3 


161 


44x16 


o 


4 


138 


44x16 


o 



thors, no effort has been made to date to quantitatively 
explain the magnitude of PMR, using the full knowledge 
of experimentally obtained sample parameters, Vb, n e , 
the mobility /i, and the quantum, or single-particle mo- 
bility /i s . Such quantitative analysis has been done in the 
present paper for ULSL samples with relatively small pe- 
riods and modulation amplitudes that allow determining 
reliable values of Vb from the CO amplitude. 20 The result 
demonstrates that magnetoresistance attributable to SO 
is much smaller than the observed PMR. We propose 
an alternative mechanism that accounts for the major 
part of PMR. After detailing the ULSL samples used in 
the present study in Sec. [H] we delineate in Sec. IIIII a 
simple analytic formula to be used to estimate the con- 
tribution of SO to PMR. Experimentally obtained PMR 
traces are presented and compared to the estimated SO- 
contribution in Sec. IIVI leading to the introduction of 
another mechanism, the contribution from drift velocity 
of incompleted cyclotron orbits, in Sec.0 which we be- 
lieve dominates the PMR for our present ULSL samples. 
Some discussion is given in Sec. lVIl followed by conclud- 
ing remarks in Sec. IVIll 



II. CHARACTERISTICS OF SAMPLES 

We examined four ULSL samples with differing pe- 
riods a, as tabulated in Table [I] The samples 
were prepared from the same GaAs/AlGaAs single- 
heterostructure 2DEG wafer with the heterointerface re- 
siding at the depth rf=90 nm from the surface, and hav- 
ing Alo.3Gao.7As spacer layer thickness of d s =40 nm. A 
grating of negative electron-beam (EB) resist placed on 
the surface introduced potential modulation at the 2DEG 
plane through strain- induced piezoelectric effect iS4 To 
maximize the effect, the direction of modulation (re- 
direction) was chosen to be along a <110> direction. 
For a fixed crystallographic direction, the amplitude of 
the strain-induced modulation is mainly determined by 
the ratio a/d. Figures ^b) and (c) display scanning elec- 
tron micrographs of the gratings. Samples 1, 3, and 4 
utilized a simple line-and-space pattern as shown in (b). 
For sample 2, we employed a patterned grating depicted 
in (c); the "line" of resist was periodically notched in ev- 
ery 575 nm by width 46 nm. The width was intended 
to be small enough (much smaller than d) so that the 
notches introduce only negligibly small modulation them- 
selves but act to partially relax the strain. The use of 
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FIG. 1: (a) Schematic drawing of the sample with voltage 
probes for measuring modulated (ULSL) and unmodulated 
(reference) part. (b),(c) Scanning electron micrographs of 
the EB-resist gratings that introduce strain-induced potential 
modulation. Darker areas correspond to the resist. A stan- 
dard line-and-space pattern (b) was utilized for samples 1, 3 
and 4. Sample 2 employed a patterned grating (c) designed 
to partially relax the strain. 



the patterned grating enabled us to attain smaller Vq 
than sample 1, which has the same period a=184 nm. 
As shown in Fig. ^a), we used Hall bars with sets of 
voltage probes that enabled us to measure the section 
with the grating (ULSL) and that without (reference) at 
the same time. Resistivity was measured by a standard 
low-frequency ac lock-in technique. Measurements were 
carried out at T—1A and 4.2 K, both bearing essentially 
the same result. We present the result for 4.2 K in the 
following. 

To investigate the behavior of PMR under various val- 
ues of sample parameters, n e was varied from about 2.0 
to 3.0xl0 15 m -2 , employing persistent photoconductiv- 
ity effect through step-by-step illumination with an in- 
frared light-emitting diode (LED). Samples 3 and 4 were 
equipped with a back gate, which was also used to alter 
n e approximately between 1.7 and 2.0xl0 15 m~ 2 . The 
electron density n e was measured by the period of CO 
or Shubnikov-de Haas (SdH) oscillation, and also by Hall 
resistivity. Concomitant with the change of n e , param- 
eters associated with the random potential scattering, fi 
and /i s , also vary. Plots of /x and /x s (the latter only for 
samples 1 and 2) versus n e are presented in Figs.|2Ib) and 
(c), respectively. Quantum mobility is deduced from 
the damping of the SdH oscillation 25 of the unmodulated 
section of the Hall bar. 

The amplitude Vq of the modulation was evaluated 
from the amplitude of CO. In a previous publication, 20 
the present authors reported that the oscillatory part of 
the magnetoresistance is given, for Vq much smaller than 



3 



0.25 



0.20 



> 

2 0.15 

E 



(a) 




J 0.02 






S(T) \J 



> 



0.00- 
120 



100 



60 

40 
30 
25- 
20 
15 



0.0 0.2 0.4 0.6 



(b) 




10- 
5 




(c) 



-+-a No. 1(184 nm) 
No. 2(184 nm) 
~ + & No. 3(161 nm) 
■^v No. 4(138 nm) 




1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 

n (10 15 m" 2 ) 



FIG. 2: (Color online) Sample parameters as a function of the 
electron density n e , varied either by LED illumination (open 
symbols) or by back-gate voltage (solid symbols), (a) Modu- 
lation amplitude Vb. (b) Mobility p. (c) Damping parameter 
pw of CO. Quantum mobility p s for samples 1 and 2 are also 
plotted by x and +, respectively. Inset in (a) shows Ap° s £ / po 
experimentally obtained by subtracting a slowly varying back- 
ground from the magnetoresistance trace (for sample 2 at 
n e =2.20xl0 15 m~ 2 , shown by solid trace) and calculated by 
Eq. @ using Vb and pw as fitting parameters (dotted trace, 
showing almost perfect overlap with the experimental trace). 




FIG. 3: (Color online) Fermi surface in the x-p x -p y space 
for (a) /3=0, (b) 0</3<l, and (c) (3=1. Each electron orbit 
is specified by the cross section of the Fermi surface by a 
constant-py plane. Streaming orbits are present in the shaded 
area. 



where A(x)=x/ sinh(x), k^J ' a = (1 / '27r 5 ! ) (akp /2)hui c with 
uj c =e\B\/m* the cyclotron angular frequency, &Q=h/e 
the flux quantum, and /Ltg=e?i/2m*(~0.864 meV/T for 
GaAs, an analogue of the Bohr magneton with the elec- 
tron mass replaced by the effective mass m*~0.067m e ). 
Apart from the factor A(tt / B) , which governs the 
damping of CO by scattering, Eq. J2J is identical to the 
formula calculated by first order perturbation theory^ 
The parameter /i\y was shown in Ref. |20| to be approx- 
imately equal to /i s , in accordance with the formula 
given for low magnetic field in the theory by Mirlin and 
Wolfle,-— Measured Ap° s £/po for the present samples are 
also described by Eq. (J2J) very well, as exemplified in the 
inset of Fig.^a). So far, we have treated the modulation 
as having a simple sinusoidal profile Vb cos(27rx/a), and 
have tacitly neglected the possible presence of higher har- 
monics. Although the Fourier transforms of Ap° s £ / p do 
reveal small fraction of the second- (and also the third- 
for samples 1 and 2) harmonics^ their smallness along 
with the power dependence on Vb of the relevant resistivi- 
ties [to be discussed later, see Eqs. (|T2*|) and (l2*2l ] justifies 
neglecting them to a good approximation. The parame- 
ters Vb and /Jw obtained by fitting Eq. J2J) to experimen- 
tal traces are plotted in Figs.|2{a) and (c), respectively. 
The latter shows piw— ^s, confirming our previous result. 
Vb does not depend very much on n e when n e is var- 
ied by LED illumination, but increases with decreasing 
n e when the back gate is used, the latter resembling a 
previous report The dependence of Vb on n e is dis- 
cussed in detail elsewhere^ Since a and d are of com- 
parable size, Vb rapidly increases with the increase of a 
(with exception, of course, of sample 2 whose amplitude 
is close to that of sample 3). Since 6<Ep <11 meV for the 
range of n e encompassed in the present study, the condi- 
tion ?7<Cl is fulfilled for all the measurements shown here 
(77=0.010-0.034). 



the Fermi Energy E F , t,=V /E f <&1, by nl - CALCULATION OF THE CONTRIBUTION 

OF STREAMING ORBITS 



Ap osc ( it \ ( T\ 

— = A I — J A \ — J In this section, we describe a simple analytic calcu- 

/3 ° \MW / \ a J lation for estimating the contribution of SO to magne- 

1 1 M Vq iff| sm [2 --^ 1 (2) toresistance. The calculation is a slight modification of 

2V27T &oPr 2 o, n^J 2 \ a J ' a theory by Matulis and Peeters^Si the theory in which 
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FIG. 4: (Color online) Functions F(/3), G{0), <j>{fi) (thin dot- 
ted, solid, and dash-dotted lines, respectively, left axis), and 
/3 2 G(/3) (thick solid line, right axis). 



semiclassical conductance was calculated for 2DEG un- 
der unidirectional magnetic field modulation with zero 
average. We modify the theory to the case for poten- 
tial modulation Vq cos(27ra:/a), and extend it to include 
a uniform magnetic field —B (the minus sign is selected 
just for convenience). The Hamiltonian describing the 
motion of electrons is given by 

s{x,p x ,Py) = [pl + {p v - eBxf] + V COS ( — 

(3) 

in the Landau gauge A=(0, — Bx, 0), where p=(j>x,Py) 
denotes canonical momentum. Using electron velocities 



de 
dp x 



Px 

m* 



de 
dp y 



Py 



eBx 



(4) 



the conductivity tensor reads (including the factor 2 for 
spin degeneracy) 



2e 2 



1 



{2irh) 2 Lx 



dx 



dp x 



df 

dp y T s ViVj(~ — ) 

V-5 Co 



(5) 

where L x represents the extent of the sample in the re- 
direction and t s denotes an appropriate scattering time 
to be discussed later >2i Since the system is periodic in re- 
direction and each SO is confined in a single period, the 
integration over x, L~ x dx, can be reduced to one 
period, a -1 dx, in calculating the conductivity from 
SO. The derivative —df/de of the Fermi distribution 
function f(e)—{l + exp[(e — EfV/cbT]} -1 ma y be a P~ 
proximated by the delta function 5(e — E-p) at low tem- 
peratures, T<S^Ep/kB- Therefore the problem boils down 
to the integration of T s ViVj over relevant part of the Fermi 
surface e(x,p x ,p y ) — E? in the (x,p x ,p y ) space. Fermi 
surface is depicted in Fig. [3] for three different values of 
f3=B/B c . Since the Hamiltonian Eq. © does not explic- 
itly include y, p y is a constant of motion that specify an 
orbit; an orbit is given by the cross section of the Fermi 



surface by a constant-p y plane. The presence of SO is 
indicated by the shaded area in Fig. |3| The ratio of SO 
to all the orbits is maximum at [3=0, decreases with in- 
creasing (3, and disappears at (3=1. 

Before continuing the calculation, we now discuss an 
ade quat e scattering time to choose. At variance with 
Ref. we adopt here unweighted single-particle scat- 
tering time T s =fi s m* /e. The choice is based on the fact 
that the angle 9=aictaxi(v x /v y ) of the direction of the 
velocity with respect to the y-axis is very small for elec- 
trons belonging to SO in our ULSL samples having small 
•q=Vo/E-p. The maximum of \0\ at a position u=2ixxja 
can be approximately written as [rjip((3, it)] 1 / 2 with 



ip(/3, u) = yl — (3 2 + (3 arcsin (3 — cos u — (3u, (6) 

whose maximum over u is given by [2770O3)] 1 / 2 with 
<j>{P)=y/l ~(3 2 +(3 arcsin (3 - (tt/2)/3, where |</>(/?)|<l for 
|/3|<1 (see Fig.0J. Since |0| is much smaller than the av- 
erage scattering angle 6> SC at~ \J 2ft s / rad estimated 
for our present 2DEG wafer, electrons are kicked out of 
SO by virtually any scattering event regardless of the 
scattering angle involved, letting r s to be the appropri- 
ate scattering time. 

The integration Eq. (J3J) over the shaded area gives the 
correction to the conductivity owing to SO, to the leading 
order in r\, as 



5a 



so 

XX 



CO 



_2_/x 
2tt 2 



/* J arcsin (3 

32V2 n. 



MP) o 

-[ W (P,u)] 3 / 2 du 



9tt 2 fi 



V 3/2 F(f3), 



(7) 



where the minus sign results because electrons trapped in 
SO cannot carry current over the (macroscopic) sample 
in x-direction and therefore should be deducted from the 
conductivity, and 



' yy 
Co 



2 p s 
2tt 2 h 

8a/2/4 

7T 2 fl 



«l(/3) 
arcsin j3 

V 1/2 G(f3), 



2[r 1 tp{(3,u)] 1/2 du 



(8) 



and Sa^.y=day2=0, where a =E F e 2 T/Trh 2 represents the 
Drude conductivity. The factor 2 in the first equalities 
accounts for the two equivalent SO areas at the upper 
and the lower bounds of p y . The functions F{(3) and 
G{(3) are defined as 



F(f3) 



16V2 



«i(/3) 



[p((3,u)fl 2 du 



arcsin /3 



and 



G(/3) = 



4%/2 



MP) 



[<p(^u)]^ 2 du, 



(9) 



(10) 



arcsin 



where the upper limit of integration ui((3) is the solu- 
tion of ip((3,ui)=0 other than arcsin (3. Both F((3) and 



5 



G((3) monotonically decrease from 1 to while /3 varies 
from to 1, as shown in Fig. 0] Since 5a^ /5ayyOC7], 
fi&xx^dvyy f° r ^^l- Correction to the resistivity by 
SO can be obtained by inverting the conductivity tensor, 
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For small 77, Sa^.® /&o can be neglected and consequently 



6pl° ^ 4V2 1 



/'B 



5/2 



T/ 5/2 

^ /3 2 G(/3). 



3/2 



(12) 



The correction therefore increases in proportion to /i, /i s , 

5/2 

and V , and decreases with a and n e . The function 
(3 2 G(f3) is also plotted in Fig. 0] which takes maximum 
at /3~0.6 and vanishes at (3=1. Our final result Eq. (fT5|) 
is identical to Eq. (41) of Ref. |2^, which is deduced for 
the case ?7<C /i s //i. (For larger 77, Ref. I23I gives somewhat 
different formula that is proportional to V 7 ^ 2 ). Note that 
our (f>((3) and G((3) are identical to the functions denoted 
as $i(/3) and & (f3), respectively, in Ref.l23t In the follow- 
ing section, Eq. (|12fl will be compared with experimental 
traces. 



IV. POSITIVE MAGNETORESISTANCE 
OBTAINED BY EXPERIMENT 

Figure [S] shows low-field magnetoresistance traces for 
samples 1—4 for various values of n e . Solid curves rep- 
resent measurements before illumination (n e varied by 
the back gate) and dotted curves are traces for n e var- 
ied by LED illumination (back gate voltage=0 V). The 
magnitude of PMR shows clear tendency of being large 
for samples having larger Vq. By contrast, the peak po- 
sitions do not vary much between samples. To facilitate 
quantitative comparison with Eq. iff^) . Fig.[5]is replotted 
in Fig. El with both horizontal and vertical axes scaled 
with appropriate parameters: the horizontal axis is nor- 
malized by B c calculated by Eq. Q using experimentally 
deduced n e and Vq shown in Fig. [2 the vertical axis is 
normalized by the prefactor in Eq. i|12|) with p s replaced 
by Pw, identifying the two parameters^ Magnetoresis- 
tance owing to SO will then be represented by a universal 
function f3 2 G((3), which is also plotted in the figures. It is 
clear from the figures that experimentally observed PMR 
is much larger than that calculated by Eq. I|12|l . Further- 
more, the peaks appear at B>B D , i.e., where SO have 




FIG. 5: (Color online) Magnetoresistance traces for various 
values of n e . Selected values of n e are noted in the figure (in 
10 15 m~ 2 ). Dotted traces indicate that the n e is attained by 
LED illumination. Note that the vertical scale is expanded 
by five times for sample 4. 



already disappeared, for all traces for samples 1—3 and 
traces with smaller n e for sample 4. The peak position 
is by no means fixed, but depends on the sample param- 
eters. This observation argues against the interpretation 
of PMR being solely originating from SO. Rather, we in- 
terpret that SO accounts for only a small fraction of the 
PMR, as suggested by Fig-EI ancl that the rest is ascribed 
to another effect to be discussed in the next section. In 
fact, humps that appear to correspond to the component 
(3 2 G(f3) can readily be recognized in traces with larger 
n e for sample 4, superposed on a slowly increasing com- 
ponent of PMR. The humps terminate at around |/3|=1, 
where the total PMR changes the sign of the curvature. 
With the increase of r]((xVo/n e ) either by decreasing n e 
(upper traces for sample 4) or by increasing Vq (samples 
1—3), (3 2 G((3) makes progressively smaller contribution 
to the total PMR, and becomes difficult to be distin- 
guished from the background. 

As has been inferred just above, the interpretation that 
the contribution Sp%® j po from SO is superimposed on an- 
other slowly increasing background component offers an 
alternative way to determine B : B e can be identified 
with the end of the hump, namely, the inflection point 
.Binf where the curvature of the total PMR changes from 
concave down, inherited from /3 2 G(/3), to concave up. To 
be more specific, B 1TL i is determined as a point where the 
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FIG. 6: (Color online) Replot of Fig. Q3 with ab- 
scissa normalized by the extinction field B c and or- 
dinate by the sample-parameter-dependent prefactor 

_____ 5/2 2 3/2 

in Eq. (1121 . apwpV a n e , with the coefficient 



a=4y2 7 r- 1 $o 1/2 MB" 5/2 =;4.04xl0 7 T 2 meV- 5/2 m- 1 . Ver- 



tical scale is expanded twice for sample 4. 
/3 2 G(/3) is also plotted for comparison. 
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FIG. 7: (Color online) (a) Magnetoresistance traces for sam- 
ple 4 with the inflection point Bi n f marked by downward open 
triangles. Traces are offset proportionally to the change in n e . 
Selected values of n e in I0 15 m -2 are noted in the figure, (b) 
Illustration of the procedure to pick up Bi n f (an example for 
n e =2.33xf0 15 m -2 ). The point at which the second deriva- 
tive (d 2 /dB 2 )(Apxx/ po) (solid curve, right axis) crosses zero 
upward (marked by open downward triangle) is identified as 
■Binf ■ Bi n f is marked also on Ap xx / Po (dotted curve, left axis). 
Shaded area indicates the contribution from SO. (c) Plot of 
Binf versus B B calculated by Eq. Q using experimentally ob- 
tained Vq. The line represents B in {=B e . 



second derivative (d 2 /dB 2 )(Ap xx / p ) changes sign from 
negative to positive as illustrated in Fig. 0b). The inflec- 
tion point i?i n f is marked by a downward open triangle 
both in (d 2 /dB 2 )(Ap xx /p ) (solid) and Ap xx /p (dot- 
ted) traces, [{d 2 / dB 2 )(Ap xx / po) shows oscillatory fea- 
tures at low field, which are attributed to the geometric 
resonance of Bragg-rcflcctcd cyclotron orbits. 33 ] Figure 
0a) illustrates the shift of B[ n { with n e . The plot of B m f 
versus B c shown in Fig. 0c) demonstrates that B in { is 
actually identifiable with B c . Thus it is now possible to 
deduce reliable values of Vb from PMR: by replacing B c 
with i?j n f in Eq. l|T)l, Unfortunately this method is ap- 
plicable only for samples with very small r/. For samples 
1—3, it is difficult to find clear inflection points because 
of the dominance of the slowly increasing component; 
(d 2 /dB 2 )(Ap xx / po) only gradually approaches zero from 
below. In the subsequent section, we discuss the origin 
of the slowly increasing background component of the 
PMR. 



V. DRIFT VELOCITY OF INCOMPLETED 
CYCLOTRON ORBITS 

An important point to be noticed is that even at the 
low magnetic- field range \B\<B e where SO is present, 
most of the electrons are in cyclotron-like orbits, namely 
the cyclotron orbits slightly modified by a weak poten- 
tial modulation, as evident in Fig. SO accounts for 
only small fraction, order of r/ 1 / 2 , of the whole orbits. 
Therefore, the contribution of these cyclotron-like orbits 
to the magnetoresistance should be taken into consid- 
eration in interpreting the PMR. We will show below 
that the slowly varying component of the PMR is at- 
tributable to the E x B drift velocity of the electrons in 
the cyclotron-like orbits that are scattered before com- 
pleting a cycle. 

It is well established that the E x B drift veloc- 
ity resulting from the gradient of the modulation po- 
tential E=— VV/(— e) and the applied magnetic field 
B=(0, 0, B) is the origin of the COii— ■ For unidirectional 
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(a) (b) 




FIG. 8: (Color online) Illustration of E x B drift velocity i> d 
affecting the electrons during the cyclotron motion. Orbits 
are depicted neglecting the modification due to the modula- 
tion V (x)=Vo cos(qx) (drifting movement and slight variation 
of the velocity depending on x) for simplicity. Top diagrams 
represent slightly larger B than bottom ones for both (a) and 
(b). On averaging v^. y along an orbit, most contribution 
comes from minimum- and maximum- a; edges, as shown by 
solid arrows in the figure. Open arrows indicate the direction 
of E x =(qVo/e)sin(qx) at the edges, (a) For B large enough 
so that electrons can travel cycles before being scattered. De- 
pending on B, va at both edges are constructive (top diagram, 
2R c /a—n +1/4 with n integer) or destructive (bottom dia- 
gram, 2R c /a—n — 1/4), resulting in maxima and minima in 
the magnetoresistance, respectively, (b) For small B so that 
electrons are scattered before completing a cycle. The inter- 
relation of Vd,y at both edges is not simply determined by B. 
The edges affect the magnetoresistance independently. 

modulation V(x)=Vq cos(qx) with q=2ir/a, the drift ve- 
locity Vd=(E x B)/£? 2 has only the y-componcnt, 

Vd, y = ^-^ sin(qx). (13) 

Electrons acquire Vd. y during the course of a cyclotron 
revolution, whose sign alternates rapidly except for when 
electrons are traveling nearly parallel to the modulation 
(#~0, 7r), i.e., around either the rightmost (maximum- a;) 
or the leftmost (minimum- a;) edges. Therefore, the con- 
tribution of the drift velocity to the conductivity comes 
almost exclusively from the two edges as depicted in Fig. 
El (a), which is actually experimentally verified in Ref. 
|35L The CO is the result of the alternating occurrence 
by sweeping the magnetic field of the constructive and 
destructive addition of the effects from the two edges, 
as illustrated by the top and the bottom cyclotron or- 
bits in Fig. [HI (a) , respectively. With the decrease of the 
magnetic field, cyclotron radius R c increases and conse- 
quently the probability of electrons being scattered be- 
fore reaching from one to the other edge increases. As a 
result, the distinction between the constructive and de- 
structive cases are blurred, letting the CO amplitude di- 
minish more rapidly than predicted by the theories^^ 



neglecting such scattering. 

The absence of CO at lower magnetic fields signifies 
that electrons are mostly scattered before traveling to 
the other edge. Although the correlation of the local drift 
velocities at the both edges is lost at such magnetic fields 
(Fig. [SI (b)) , each edge can independently contribute to 
the conductivity. It is to this effect that we ascribe the 
major part of PMR in our ULSL samples. Note that the 
onset of CO basically coincide with the end of the PMR, 
bolstering this interpretation. 

It can be shown, by an approximate analytic treatment 
of the Boltzmann's equation, that the effect actually gives 
rise to PMR with right order of magnitude to explain the 
experimentally observed slowly-varying component. For 
this purpose, we make use of Chambers' forrnulaj2£i2LiS 
representing the relaxation time approximation of Boltz- 
mann's equation, to obtain, from the drift velocity, the 
component D yy of the diffusion tensor, 

/>oo 

Dyy = / e^(v^ y (t)v d .jO))dt, (14) 
Jo 

where (...) signifies averaging over all possible initial con- 
ditions for the motion of electrons along the trajectories. 
Einstein's relation is then used to obtain corresponding 
increment in the conductivity, da y y=e 2 D(Ep)Dyy with 
D(EF)—m*/'Kh 2 =(^ol J -B) 1 the density of states, and fi- 
nally it is translated to the resistivity by tensor inversion, 
Spxx/po—i^cT^Sayy/ao- We use unperturbed cyclotron 
trajectory, x=X + R c cos 9, for simplicity, neglecting the 
modification of the orbit by the modulation (and accord- 
ingly, SO is neglected in this treatment), which is justified 
for small rj. Since the initial condition can be specified 
by the guiding center position X and the initial angle 9q, 
we can write 

K,(tK,(o))=(^) 2 ^^^/; fo 

sin{q[X + R c cos(6» + uj c t)}} sin[q(X + R c cos0 o )].(15) 

Therefore Eq. (|14f) can be rewritten, performing the in- 
tegration over t first, as 

Dyy= (S) lJ Q dX ^ J rf 0osm[g(X+i?cCos0o)]/(0o), 

(16) 

with 

I(0 O ) = / e ~ t/T sm{q[X + R c cos(6» + uJ c t)]}dt. (17) 
Jo 

Evaluation of Eq. I|16|) for a large enough magnetic field 
reproduces basic features of Eq. (j2J , as will be shown in 
the Appendix. Here, we proceed with an approximation 
for small magnetic fields. The approximation is rather 
crude but is sufficient for the purpose of getting a rough 
estimate of the order of magnitude. 

Because of the exponential factor, only the time t<r 
contributes to the integration of Eq. (|17|) . Due to the 
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rapidly oscillating nature of the sin{} factor and the 
smallness of oj c t<oj c T, I(9 Q ) takes a significant value only 
when 9q resides in a narrow range slightly below — or 
~7r, corresponding to the situation when electrons trav- 
els near the rightmost or the leftmost edge, respectively, 
within the scattering time. It turns out, by comparing 
with the numerical evaluation of Eq. (|17fl using sample 
parameters for our present ULSL's, that the following 
approximate expressions roughly reproduce the right or- 
der of magnitude and the right oscillatory characteristics 
(the period and phase) of Eq. H17fl for low magnetic field 
(|J3|<0.02 T): 



Wo) 



TTr[sin(qX)J (qR c 
TTr[sm(qX)J (qR c 



cos(qX)H (qR c )} (# ~0) 
co S (qX)H (qR c )} (8 a ~ir) 

where Jq(x) and Ho (a;) represent 0-th order Bessel and 
Struve functions of the first kind, respectively. The ap- 
proximation can be obtained by replacing the exponential 
factor by a constant lo c t and limiting the range of the 
time integral to include only one edge. Here we noted 
that the integration of cos(qR c cos9) and sin(gi? c cos 9) 
over the range of 9 including either of the rightmost (9=0) 
or the leftmost (9=n) edge can be approximated (since 
only the close vicinity of the edges makes significant con- 
tribution to the integration) by, 




o.io 



d9 

rightmost 



r/2 



d9, 



I-tt/2 

and used the relations 

r Tr/2 



d9 

leftmost 



3tt/2 



d9, (19) 



FIG. 9: (Color online) Replot of magnetoresistance traces 
normalized by the prefactor in Eq. <Y2'2l . <■ 



tt/2 



a' n V a "n e 

with a' = (7r/2) 1/2 $o Vb" 2 =4.06x10 14 T meV -2 m" 2 , after 
subtracting the contribution from SO, 8p^/po in Eq. 1121 . 
Contribution attributable to drift velocity of incompleted cy- 
clotron orbits is given by \B\, which is also plotted by dash- 
dotted line. 



cos(qR c cos 9) d9 = / cos(qR c cos 9)d9 = tt Jo(qR c 

-tt/2 J-k/2 



and 

k/2 

sm(qR c cos9)d9 

n/2 



3tt/2 

sin(qi? c cos 

/2 



TrH Q {qR c ). 

(20) 



Substituting Eq. JTSJ) to Eq. (UHJ results in 



D 



and with Einstein's relation one finally obtains 



g drift 



Pa 



1 fi 2 V 2 



2 $oA*£' 



3/2 



\B\ 



(21) 



(22) 



Here we made use of asymptotic expressions Jo(x) w 
{2/nx) 1 / 2 cos(x-tt/4) and H (a;)w(2/ irx) 1 / 2 sin(a; — 7r/4) 
valid for large enough x (corresponding to small enough 
B). 

In order to compare experimentally obtained PMR 
with Eq. l(2^|) , PMR traces shown in Fig. are replot- 
ted in Fig. [5] normalized by the prefactor in Eq. (1221) , 
after subtracting the small contribution from SO repre- 
sented by Eq. I)12p. The scaled traces show reasonable 



agreement with \B\ at low magnetic fields, as predicted 
in Eq. I)22[l. testifying that the mechanism considered 
here, the drift velocity from incompleted cyclotron or- 
bits, generates PMR having the magnitude sufficient to 
explain the major part of PMR observed in our present 
ULSL samples. Possible sources of the remnant devi- 
ation, apart from the crudeness of the approximation, 
are (i) the neglect of higher harmonics and (ii) the ne- 
glect of negative magnetoresistance (NMR) component 
innate to GaAs/AlGaAs 2DEG 39 arising from electron 
interactions^Si^i or from semiclassical effecti^Sil 3 . The n- 
th harmonic gives rise to additional contribution analo- 
gous to Eq. with Vq and a replaced by the amplitude 
V n of the n-th harmonic potential and a/n, respectively, 
and therefore, in principle, enhances the deviation. In 
practice, however, the effect will be small because of the 
small values of V n and its square dependence. On the 
other hand, the discrepancy can be made smaller by cor- 
recting for the NMR. We have actually observed NMR, 
which depends on n e and temperature, in the simultane- 
ously measured "reference" plain 2DEG adjacent to the 
ULSL (see Fig.[T](a)). Assuming that the NMR with the 
same magnitude are also present in the ULSL part and 
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superposed on the PMR (the assumption whose validity 
remains uncertain at present), the correction are seen to 
appreciably reduce the discrepancy. 

The approximation leading to Eq. 12211 is valid only 
for very small magnetic fields. With the increase of 
the magnetic field, the cooperation between the leftmost 
and the rightmost edges is rekindled, and the magne- 
toresistance tends to the expression appropriate for large 
enough magnetic field, outlined by Eq. (|A.11|) . which 
includes a non-oscillatory term (the first term) as well 
as the term representing CO (the second term). Note 
that the non-oscillatory term approaches a constant, 
a'yU 2 V^ 2 a -1 7j,r 3 / 2 (m*/27reT), at small magnetic field, al- 
though the exact value of the constant is rather diffi- 
cult to estimate due to the subtlety in choosing the right 
scattering time r, as will be discussed in the Appendix. 
Therefore the (linear) increase of Sp x " lt /po with \B\ is 
expected to flatten out at a certain magnetic field. The 
peak in the PMR roughly marks the position of this tran- 
sition, which basically corresponds to the onset of the 
cooperation between the two edges. Thus the peak po- 
sition is mainly determined by the scattering parame- 
ters and is expected to be insensitive to Vb, in agreement 
with what has been observed in Fig. [3] Experimentally, 
the peak position B p is found to be well described by 
an empirical formula B p (T)=[4-\/2^w(ni 2 /Vs)] _:L , using 
^t\y determined from CO. On the other hand, the height 
of the PMR peak are seen to roughly scale as Vb 2 , as 
inferred from Fig. [5] which reveals that the normalized 
peak height tends to fall into roughly the same value (no- 
tably the top panel showing two samples having the same 
a and different Vq), so long as the period a are the same. 
This is better shown after correcting for the NMR effect 
mentioned above. The height of the normalized peak 
slightly decreases with decreasing a (roughly proportion- 
ally to a), resulting in an empirical formula for the peak 
height (Ap xx /p )v cak ~ 3xl0- 3 [^(m 2 /Vs)] 2 [U (meV)] 2 
[n e (10 15 m~ 2 )]- 3 / 2 . (Unfortunately, sample 4 with larger 
n e significantly deviates from this formula.) 



VI. DISSCUSSION ON THE RELATIVE 
IMPORTANCE OF THE STREAMING ORBIT 

Although PMR was thus far generally interpreted to 
originate from SO, contribution from mechanisms other 
than SO was also implied in theoretical papers. By 
solving Boltzmann's equation numerically, Menne and 
Gerhardts 21 calculated PMR and showed separately the 
contribution of SO which did not account for the entire 
PMR (see Fig. 4 in Ref. l2l|) . leaving the rest to alter- 
native mechanisms (although the authors did not discuss 
the origin futher) . Mirlin et al£& actually calculated con- 
tribution of drifting orbit, which is basically similar to 
what we have considered in the present paper. They pre- 
dicted cusp-like shape for the magnetoresistance arising 
from this mechanism, which is not observed in experi- 
mental traces. In both papers, the major part of PMR is 



still ascribed to SO, with other mechanisms playing only 
minor roles. In the present paper, we have shown that 
the relative importance is the other way around in our 
ULSL samples. However, we would like to point out that 
the dominant mechanism may change with the amplitude 
of modulation in ULSL. 

The reason for the contribution of SO being small in 
our samples can be traced back to the small amplitude 
of the modulation, combined with the small-angle na- 
ture of the scattering in the GaAs/AlGaAs 2DEG. As 
mentioned earlier, small r]=Vo/E-p limits the SO within 
narrow angle range \0\<y/2rj, letting the electrons being 
scattered out of the SO even by an small-angle scatter- 
ing event, hence the use of r s in Eq. This leads to 
small S<7yy , since t s <^t. Within the present framework, 
relative weight of SO in PMR decreases with increasing 
77, since the ratio of Eq. (|12(l to Eq. (|22|l is proportional 
to ?7 _1 ' 2 , in agreement with what was observed in Fig. 
However, the situation will be considerably altered 
with further increase in r\ (typically 7?>0.1). Then, due 
to the expansion of the angle range encompassed by SO, 
electrons begin to be allowed to stay within SO after 
small-angle scattering, requiring r s in Eq. JSJ to be re- 
placed by larger (possibly -B-dependent) values. In the 
limit that the range of \Q\ is much larger than the av- 
erage scattering angle, r s should be supplanted by ordi- 
nary transport lifetime (momentum-relaxation time) r, 
resulting in much larger 5(Tyy . This largely enhances the 
relative importance of SO, possibly to an extent to ex- 
ceed the contribution from the drift velocity. We presume 
that the contribution of SO is much larger than in our 
case in most of the experiments reported so far which 
showed the shift of PMR peak position with the modula- 
tion amplitude^ii^ii^iSiii. Even in such situation, how- 
ever, it will not be easy to obtain simple relation between 
the peak position B p and the amplitude Vb because of the 
complication by the remnant contribution from the drift 
velocity. In most experiments, Vb is varied by the gate 
bias, which concomitantly alters the electron density and 
scattering parameters, thereby affecting the both contri- 
butions as well. 



VII. CONCLUSIONS 

The positive magnetoresistance (PMR) in unidirec- 
tional lateral superlattice (ULSL) possesses two different 
types of mechanisms as its origin: the streaming orbit 
(SO) and the drift velocity of incompleted cyclotron or- 
bit. Although virtually only the former mechanism has 
hitherto been taken into consideration, we have shown 
that the latter mechanism account for the main part of 
PMR observed in our ULSL samples characterized by 
their small modulation amplitude. The share under- 
taken by SO decreases with increasing r)=Vo/E-p, inso- 
far as r\ is kept small enough for the electrons in SO to 
be driven out even by a small-angle scattering character- 
istic of GaAs/AlGaAs 2DEG; 7y<0.034 for our samples 
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fulfills this requirement. In this small rj regime, the peak 
position of PMR is not related to the modulation ampli- 
tude Vb but rather determined by scattering parameters; 
the peak roughly coincide with the onset of commensu- 
rability oscillation (CO) that notifies the beginning of 
the cooperation between the leftmost and the rightmost 
edges in a cyclotron revolution. The height of the peak, 
on the other hand, are found to be roughly proportional 
to Vq. For small enough 77, the contribution of SO be- 
comes distinguishable as a hump superposed on slowly- 
increasing component and the magnetic field that marks 
the end of the SO, B c , can be identified as an inflection 
point of the magnetoresistance trace where the curvature 
changes from concave down to concave up. The extinc- 
tion field B e provides an alternative method via Eq. (JJJ 
to accurately determine Vq. We have also argued that for 
samples with 77 much larger than ours, typically 77^0. 1, 
the relative importance of the two mechanisms can be 
reversed and the PMR peak position B p can depend on 
Vo, although it will be difficult to deduce a reliable value 
of Vb from B p . 
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APPENDIX: APPROXIMATION FOR HIGHER 
MAGNETIC FIELD 

In this appendix, we delineate the approximation of 
Eq. (|14|l at higher magnetic field, which leads to an ex- 
pression for commensurability oscillation (CO). When 
the velocity Ud,j/(t) is a periodic function of time with 
period T, Eq. (jfij reduces to2& 

1 f T 

D vv = 1 _ e _ T/T / e~ t/r < Vd,y(*)»d,w(0) > dt. 

(Al) 

Using here again the unperturbed cyclotron orbit x= X + 
R c cos(w,4), one obtains 

Dyy = (SI) lj dx h J de o Mq(x + R c cos e )) 

1 _l-T/r I T(°o) 

(A.2) 

with T=2tt/lo c and 

r T 

I T (0 O ) = / e ~ t/T sm{q[X + R c cos(6 +uj c t)]}dt. (A.3) 
Jo 



Again because of the sin{} factor, the main contribu- 
tion in the integration comes from the narrow band of t 
around 9$ + uj c t^0 (or 2tt depending on the initial angle 
9o) and tt. For large enough lu c , the band of t becomes 
narrow enough to allow the exponential factor e~ l l T to 
be approximated by a constant value at t=—(&o — kir)/u> c 
(with fc=0,l, and 2). Thus, using the relations Eq. J20J, 
It(@o) can be approximately written, depending on the 
values of 9q, as 

TT e 

I<(9 ) ^ — e— [sm(qX)J (qR c ) + cos(qX)U (qR c )} 

UJ C 

+— e —[sm(qX)Jo(qR c ) - cos{qX)U (qR c )} 

(A.4) 

for -7T + <# <0~ and 

I>(9 ) —e^[sm(qX)J (qR c ) - cos(qX)H {qR c )] 

+—e 2 ^ L [sin(qX)J {qR c ) + cos(qX)n a (qR c )} 

(A.5) 

for + <9o<n~ , where the superscripts + and — indicate 
small setbacks to avoid the region where the integration 
have significant value. When #0 approaches closer to the 
boundary, It(0q) approaches the average of the values on 
both sides, 

It(0o -> 0) - [I<(9 )+I>(9 Q )}/2 (A.6) 

and 

I T (9 -^t:)^ [I>(0 Q ) + I<(e o - 2tt)]/2, (A.7) 
which can be shown by using the relations, 

r° r' 2 tt 

I cos(qR c cos 9)d9 = I cos(qR c cos 9)d9 = — J (qR c ) 
J-tt/2 Jo 2 

(A.8) 

and other related equations corresponding to the halves 
of Eq. . In the integration by 9 in Eq. (|A.2(I , only 
^0^0 and tt contributes to the integral for the same rea- 
son as before. The integration, after slightly shifting the 

limits of the integral from to + , yields 

TT 

{(1 + e~&) [sm 2 (qX)J 2 (qR c ) + cos 2 («X)H 2 (qR c )] 

+2e-^ [sin 2 (qX), J 2 (qR c ) - cos 2 (qA)H 2 (gi? c )] } . 

(A.9) 

Finally, by averaging over X, Eq. (|A.2(1 becomes 
vv 2\eB J 2w c 
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jcoth [Jo 2 {qR c ) + H 2 (qRc)] 

+ sinh -i fjL) [j 2 (qR c )~H 2 (qR c )] 



eB 



1 



• coth 



irqR c 
n 

LO c T 



A 



LU C T 



sin(2qR c 



(A.10) 



which can be translated to magnetoresistance with the 
use of Einstein's relation, resulting in, 



1 



1 



Pa 



2V2tt $oMb' 
— coth 



a 

IT 
UJ C T 



3/2 



\B\ 



A 



7T 

UJ C T 



sm(2qR c 



(A.ll) 



The formula agree with the asymptotic expression of Eq. 
(21) in Ref. |23 for large enough magnetic fields. The 
second term in Eq. (|A.11|) represents CO, which repro- 
duces qualitative features of Eq. ((5J). It should be noted, 
however, that the anisotropic nature of the scattering in 



GaAs/AlGaAs 2DEG cannot be correctly treated in the 
present simple relaxation-time-approximation approach 
employing only one scattering time. The anisotropic 
scattering plays an important role in CO because of its 
high sensitivity to the small-angle scattering. There- 
fore Eq. (|A. 1 1|) is of limited validity to describe CO 
in our ULSL. However, it is interesting to point out 
that if we are allowed to replace r only in the scat- 
tering damping factor A(tt/lu c t) with the single-particle 
T s =/i s m*/e~/xw'Ti*/e, we acquire Eq. (j2J except for the 
thermal damping factor. The choice of r s for the damping 
is not unreasonable, considering that the factor A(tt/lu c t) 
stems from the exponential factor in Eq. ljA.9l) . which 
describe the cooperativeness between the leftmost and 
rightmost edges that is susceptible to a small-angle scat- 
tering. [In general, velocity-velocity correlation in Eq. 
for t^kT/2 with fe=l,2,3,..., namely between the 
edges separated by more than half revolution of the cy- 
clotron orbit, requires precise positioning after the revo- 
lution, which is ruined by small angle scattering. There- 
fore the use of t s is reasonable. However, this does not 
justify the replacement only in the damping factor.] The 
thermal damping factor A(T/T a ) can readily be incor- 
porated by allowing for thermal smearing of the Fermi 
edge, namely by taking the average over the energy of 
the sm(2qR c ) term weighted by the factor (—df/de). 
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